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The thermodynamic instabilities of a binary mixture of sticky hard spheres (SHS) in 
the modified Mean Spherical Approximation (mMSA) and the Percus-Yevick (PY) 
approximation are investigated using an approach devised by X. S. Chen and F. 
Forstmann [J. Chem. Phys. 97, 3696 (1992)]. This scheme hinges on a diag- 
onalization of the matrix of second functional derivatives of the grand canonical 
potential with respect to the particle density fluctuations. The zeroes of the small- 
est eigenvalue and the direction of the relative eigenvector characterize the instability 
uniquely. We explicitly compute three different classes of examples. For a symmet- 
rical binary mixture, analytical calculations, both for mMSA and for PY, predict 
that when the strength of adhesiveness between like particles is smaller than the one 
between unlike particles, only a pure condensation spinodal exists; in the opposite 
regime, a pure demixing spinodal appears at high densities. We then compare the 
mMSA and PY results for a mixture where like particles interact as hard spheres 
(HS) and unlike particles as SHS, and for a mixture of HS in a SHS fluid. In these 
cases, even though the mMSA and PY spinodals are quantitatively and qualitatively 
very different from each other, we prove that they have the same kind of instabil- 
ities. Finally, we study the mMSA solution for five different mixtures obtained by 
setting the stickiness parameters equal to five different functions of the hard sphere 
diameters. We find that four of the five mixtures exhibit very different type of insta- 
bilities. Our results are expected to provide a further step toward a more thoughtful 
application of SHS models to colloidal fluids. 
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I. INTRODUCTION 



Thermodynamic instabilities are important to locate on the phase diagram of a fluid 
system those regions where the system can not exist as a single phase. 

For a one-component system with Helmholtz free energy A, pressure P, in a volume V, 
at a temperature T, the condition for phase stability is (d 2 A/dV 2 ) TN = — (dP/dV) T N = 
1/ (Yxt) > 0. The points where the isothermal compressibility \t diverges define the so 
called spinodal line, or phase instability boundary that separates the stable from the 
unstable region of the phase diagram. In the stable region, where \t > 0, the system can 
exist in a single phase, while inside the other region the free energy can be lowered by 
phase separation into two phases with different densities. This kind of instability is usually 
referred to as mechanical instability, associated with a gas-liquid transition or condensation 

In a binary mixture the situation is more complex jjj, |2j, y, |4| , y, |y] . The thermodynamic in- 
stability is located on the points of the phase diagram where (d 2 G/dx 2 ) T p N /xt = 0, where 
x is the concentration of one of the two species, and G is the Gibbs free energy. The points 
where Xt* 1 = are instabilities of pure condensation (and the Bhatia Thornton Q density- 
density structure factor, S pp {k) , diverges at k = 0). The points where (d 2 G/dx 2 )T,p,N = 
are again instabilities of pure condensation when 6 = p(t>i — v 2 ) = (dV/dx) T p N /V diverges 
(p is the total number density, Vi the partial molar volume, per particle, of species i. In 
this case all Bhatia Thornton structure factors diverge at k — 0) and are instabilities of 
pure demixing when 5 = (in this case the Bhatia Thornton concentration-concentration 
structure factor, S xx (k) diverges at k — 0). But, in general (for an asymmetric mixture), the 
kind of instability may be in between one of pure condensation and one of pure demixing, 
with 5 finite and different from zero (also in this case all Bhatia Thornton structure factors 
diverge at k — 0). For the particular case of a binary symmetric mixture the only allowed 
instabilities are the ones of pure condensation and of pure demixing, since 6 = 0. 

A different route was followed by Chen and Forstmann p to characterize the instability 
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uniquely in terms of an angle a, function of the density and x. 

The purpose of this work is to investigate the nature of instabilities for a binary mixture 
of sticky hard spheres (SHS). The SHS one-component model was originally proposed by 
Baxter j^, 0, 3], who showed how it admitted an analytic solution in the Percus-Yevick 
(PY) approximation. The PY solution was later extended to mixtures 
and it is nowaday regarded as extremely useful in colloidal systems. In the SHS model one 
accounts for a very short range attractive potential by defining an infinitely narrow and deep 
square well. This limit is carried out in a suitable way so that the second virial coefficient 
is finite. Due to its highly idealized nature, the one-component SHS model is not free of 
pathologies 1151. Nonetheless this model has recently regained considerable attention in 

nnnn 

studies of colloidal suspensions [la, 1121 LUa, H2J especially in its polydisperse version. Since 
the PY solution of a p— component SHS mixture requires the solution of p(p + l)/2 coupled 
quadratic equations which are hard to solve for high p, attempts have been made to treat 



the model with "simpler" approximations 



2 it ] . which would allow analytic solution even 



for polydisperse systems. One of these approximations, that we will consider in this work, 
is the modified Mean Spherical Approximation (mMSA) 22]. 

In the present work we apply the Chen and Forstmann formalism to a binary SHS mixture, 
using both the mMSA and the PY approximation. The former can be regarded as the zero 
density limit of the latter and hence its predictions must be accepted with care. However 
it has its main merit in the fact that it entitles analytical predictions even in the multi- 
component case, unlike the PY closure. 

Three classes of systems will be discussed in details. First we consider the symmetric 
mixture, where equal-size equimolar components interact with variable strength only in the 
unlike part. This simplified case was already studied by Chen and Forstmann for hard 
core particles with attractive Yukawa interactions within the reference hypernetted chain 
approximation. In this particular SHS case we are able to perform a full characterization of 
the mixture both for mMSA and PY. In a second class we discuss two paradigmatic cases: 
(i) A fluid having HS interactions among like particles and SHS interactions for the unlike 
(System A) and (ii) A fluid formed by one SHS species and another HS one (System B). For 
PY both cases have been previously discussed by Barboy and Tenne ^J], by Penders and 
Vrij [23], and by Regnaut, Amokrane, and Heno 0,13 without, however, tackling the issue 
of the stability nature. Even in these two detailed analytical investigation can be 
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carried out. Building upon our recent work 



26| , we finally discuss a third class of examples 



involving a general binary mixture where, however, the stickiness parameters are related to 
the sizes of the particles according to some plausible prescriptions . Within the mMSA, 
we are then able to discuss the nature of the instabilities previously calculated in Ref. 
by evaluating numerically the Chen and Forstmann angle a. 

The remaining of the paper is organized as follows. In Section |H] we briefly outline Chen 
and Forstmann's approach, in Section IIHI we report the PY and mMSA solutions for the 
Baxter factor correlation function of the SHS mixture. Section ITVl is dedicated to the binary 
symmetric mixture, whereas Section|3to Systems A and B. Section IVll deals with five binary 
mixtures obtained setting the stickiness parameters equal to five different functions of the 
sphere diameters. 



II. METHOD FOR ANALYZING THE INSTABILITY 



For the sake of completeness, we briefly recall the main steps of the method reported in 
Ref. [5]. In doing this, however, we shall follow the general density functional formalism 
outlined in Ref. |27[ which yields a clearer viewpoint. 



A. The Chen and Forstmann formalism 

Consider a binary mixture with N\ particles of species 1 with coordinates r\,...,r 1 Nl 
and N2 particles of species 2 with coordinates rf, . . . , r 2 N2 interacting through spherically 
symmetric pair potentials. Define the microscopic densities to be 

p,(r) = ^(r-rt) t = 1, 2 (1) 

for each one of the two species. 

Consider now the non- homogeneous system with an external potential 0i(r) acting on 
the particles of species 1 and an external potential 02 ( r ) acting on the particles of species 2. 
Let Hi and Aj be the chemical potential and the de Broglie thermal wavelength, respectively, 
for species i, N = N\ + N 2 the total number of particles, and = ({r^}, {r^}) a short-hand 
notation for the total set of coordinates. The grand partition function of the system with 
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total internal energy W(r N ) is a functional of the generalized potentials tij(r) = (3[pi — 0j(r)] 

OO OO If 

Ar n n A? A^l! An 2 ]V 2 ! J 



e -/3n[ui,« 2 ] ^2) 



where f2 is the grand free energy. It can be proven |27( that the functional Q is strictly 
concave in u\ and u 2 (if we opportunely restrict its domain of definition). The equilibrium 
number density of species % is given by 

Pi (r) = fa r = . 3 

It follows that the following functional of {p{\ and {u{\ 

pA[p 1 ,p 2 ,u 1 ,u 2 ] = ^ y pi(r)ui(r) rfr + /3ft [ui, w 2 ] , (4) 



»=i 



is also strictly concave in wi and u 2 , so it admits a unique maximum for it, — Ui, i — 1,2, 
where the {Mj} can be determined univocally from Eq. (J3J) once the equilibrium densities 
are known. 

- n 

We now set A[p 1; p 2 ] = A[pi,p 2 , Ui,u 2 }. Again one can prove |27( that this Helmholtz free 
energy is a strictly convex functional in p\ and p 2 . 

Introduce the following "grand free energy functional" of the densities 

2 r 

(3Q'[ Pl , p 2 \ = (3A[ Pl , P2 ]-J2 / P*(rK(r) dr , (5) 

i=i J 

where {vi} are some given generalized potentials, independent of the densities. Clearly only 

when Vi = iti, i = 1, 2, we have Q' = Q, i.e. equilibrium. 

Taking the first functional derivative of Q' with respect to the densities we find 

PWj SpA[ Pl ,p 2 ] 

5 Pl {v) 5 Pl {v) V ' l[r) 

= Ui(r) - Vi(r) , (6) 

where in the second equality Eqs. (jHJ) and where used. At equilibrium we then have that 
the first functional derivatives of Q' vanish and ft' attains its minimum value. 

The second functional derivatives of ft' with respect to the densities at equilibrium are 



W[P1,P2] 



_ fatj(ri) 

equil. 5 Pi^) 



_ SjjSfa - r 2 ) 

~~ 777^1 Q^r^raj , {() 

equil. PA L l) 
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where Cjj(ri,r2) are the partial direct correlation functions of the system. 

So a Taylor expansion, up to the second order terms, yields the fluctuation of fl' around 
the equilibrium caused by small density fluctuations 



1 




E 



5 ij 5{r 1 - r 2 ) 



- c ij (r 1 ,r 2 ) 



SpiirijSpjfa) drxdr-z . 



(8) 



If the system is homogeneous and isotropic at equilibrium (i.e. «j(r) = (3pi, i = 1,2), so 
that 



V 



Pi 



c ij (r 1 ,r 2 ) = CiA\ri -r 2 |) 



(9) 
(10) 



where V is the volume (assumed large enough), then we can rewrite the integral of Eq. (jHJ), 
which is a convolution, as a k-integral of a product of Fourier transforms. Replacing the 
k-integral [(27r)~ 3 J dk. . .1 by a sum over discrete k-values [V^ 1 . . .], one obtains 

1 1 



(ii) 



where Spi(k) = 5pi(k)/ \fpl and the asterisk indicates complex conjugation, having denoted 
with the tilde the Fourier transform 



so that 



/(k) = / /(r)e^ r dv , 
>v 



(12) 



(13) 



Notice that, due to the symmetry of the direct correlation functions under exchange of 
species indexes, the matrix A(k) is symmetric. 

The probability distribution for the density fluctuations 5pi (at constant T, V, and {pi}) 
is proportional to e~^ sn ' 
products 



28| . We therefore get for the mean values of the fluctuation 



<*tf(k)tfft(k)> = V 



-i 



(k) 



V[Sij + y/plpj hij(k)] 



(14) 
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where the last equality exploits the Ornstein-Zernike (OZ) equations between the partial 
total correlation functions hij and the partial direct correlation functions. 

Next define the molar fraction of species i to be Xi = Pi/p, with p = J2iPi being the 
total density of the mixture. One usually introduces |7j two linear combinations of fluctua- 
tions of partial densities, i.e. the fluctuation of total density, Sp(k), and the fluctuation of 
concentration of species 1, 8x(k), 

5p(k) = 6p 1 (k) + 6p 2 (k) (15) 

= \M\/^Pi( k ) + v / ^2^P 2 (k)] , 

8x{k) = \[ P2 5~ Pl (k) - Pl 5p 2 (k)} (16) 
P 



^^-[ y /x 2 '5p 1 (k) - ^/x[5p 2 (k)} , 

so that, if 8p\ and 5p 2 change in proportion to their respective mean concentration, then 
5x = 0. 

We also introduce 

5p(k) = ^P(k) , (17) 
Sx(k) = x [^-5x(k) , (18) 

V X1X2 

so that, in terms of the following two column vectors 

u^M, ^-( Sm V (19) 
\5p 2 (k)J \5x(k)J 

Eqs. (fT5^) and (fTH|) can be written in compact notation as u = Uv where 



U=| ^ ^ I , (20) 



/X 2 -yJX-i 

notice that U 2 = I, where I is the identity matrix. 

We find then from Eq. (|11|) (superscript T indicating the transpose) 

M = ~^v T *(k)M(k)v(k) , (21) 
where M(fc) is the following symmetric matrix 

M(fc) = UA(A;)U =| PP PX I , (22) 

M xp M xx 



8 



with 

M pp = 1 - p[x\c n + x\c 22 + 2x 1 x 2 c 12 ] , (23) 

M xx = 1 - pxix 2 [cn + c 22 - 2ci 2 ] , (24) 

M px = = Pv/x^X2"[x 2 C22 - xicn - (x 2 - xi)ci 2 ] . (25) 

The elements of the M(0) matrix are related to thermodynamic quantities as shown in 
the Appendix. In particular the determinant of M is 

det(M )^ 2 <^gf)^, (26) 

where xt is the isothermal compressibility, and Xt = PI P ls ^ ne isothermal compressibility 
of the ideal gas. 

For the particular systems that we shall consider in the following, it turns out that the 
matrix A can be written, using the Wiener- Hopf factorization in terms of the Baxter factor 
matrix Q jl]] 

A(k) = Q T *(k)Q(k) . (27) 

Hence det[M(fc)] = det[A(fc)] = | det[Q(A;)]| 2 > and trace[M(A;)] = trace[A(Jfe)] > 0. 

The inverse of M(A;) yields the mean square fluctuations of total density and concen- 
tration, i.e. the density- density structure factor Sn P {k), the concentration-concentration 
structure factor S xx (k), and the cross term S px (k) jj| 

S PP (k) = i<^(k)tfp(k)> = [M-\ p (k) , (28) 
S xx (k) = ^<5x*(k)(J5(k)> = x,x 2 [M-% x (k) , (29) 

S P .(k) = y^(6p*(k)6x(k)} = ^T 2 [M-\ x {k) . (30) 

Now, since M(fc) is a symmetric matrix, it can be diagonalized through an orthogonal 
change of basis and it will have real eigenvalues 

„. trace [M(fc)] ± J {tmce\M(k)}} 2 - 4 det[M(fc)] , . 

^±{k) = ^ ' ( 31 ) 

with \+(k) > \~{k) > 0. For the normalized eigenvectors we find 

z ± (k) = ( a±{k) | , (32) 
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with 

The transition matrix to the base formed by the eigenvectors will be 

Z(k) = I ^ a ~ {k) I , (35) 



b+(k) b-(k) 



Eq. f!21l) can then be recast into the form 

1 1 

2pV 



Sn'(5 Pl ,8p 2 ) = ~Y,lK(k) |5p + (k)| 2 + A_(*) |5p_(k)| 2 ] , (36) 



where 5p± are the Fourier components of the vector for the total density and concentration 
fluctuation in the eigenvector base, namely 

Z-'v = | S?+ I , (37) 



or 



5p_ 



Sp + (k) = a + {k) Sp{k)+b + (k) Sx (k) , (38) 
5p_(k) = a_(Jfc) 5p(k) +6_(A;) 5x (k) . (39) 



B. Characterization of the instability 

We wish to know which combination of density and concentration fluctuations, (5p, Sx) 
or (5p + , 5p~-), yields the smallest increase 5ft 1 of grand free energy. The border of a stability 
region (spinodal line) will be determined by the smaller eigenvalue X-(k) going to zero. It 
is important to remark that the minimum eigenvalue will vanish if and only if det[M(/c)] = 
X-(k)X + (k) = | det[Q(A;)]| 2 vanishes. The spinodal equation thus corresponds to 

A_(Jfe) = or det[Q(fc)] = 0. (40) 

For all k- vectors with k = |k| being a solution of the spinodal equation, we can cal- 
culate the related eigenvector z_(fc) and find, from Eq. (|39|). one non-zero linear com- 
bination <5p_(k) of density and concentration fluctuations for which 5Q' = 0. Thus 
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n = ai'ctai) ( — j = arctan 



(41) 



z_(fc) = [a_ (A;), £>_(&)] characterizes the phase transition uniquely. On defining the angle 
(see Fig. P) 

Q 12(k) - Qll(k)y/X2 
= arc l an — — — — — — 

k=k \_Ql2{h)y/x^+Qii(k)^/x{ 

the instability will be predominantly of the demixing type when a is close to (i.e. only 
concentration fluctuations occur) and predominantly of the condensation type when a is 
close to ±7r/2 (i.e. only density fluctuates at fixed concentration). 

The same feature can be seen in real space. When X-(k) = and a = (=> a_ = 0, b + = 
0, and therefore 5p + = a + 8p, 5p- = 6_ 5x), one can get 5Q' = only if 5p+(k) = 0, which 
requires Sp(k) = 0, i.e. the fluctuations that do not increase the "grand free energy" can be 
expressed as 

k 

|k]=fc 

5p 2 (r) = -S Pl (r) . (43) 

On the other hand, when \-(k) = and a = ±7r/2 (=>• 5p + = b + Sx, 5p- = a_ Sp), 
5p + (k) = now requires &r(k) = 0, which corresponds to 

5p 2 (v) = +^5 Pl (v) . (44) 
Pi 

Eq. (J4*2*|) yields oscillating partial density fluctuations for species 1 on the spinodal line, 
whereas Eqs. ()43j) and (|44|) represent the two different behaviors of the species 2 in cor- 
respondence to the two extreme values of a (0 and ±7r/2, respectively). For a = 0, the 
fluctuations of species 2 must be in opposition of phase compared to those of species 1, 
[see Eq. (JIEJ)], and this can be clearly interpreted as related to spatial demixing. In the 
opposite case (a = ±7r/2), Eq. (jUj) means that an increase of species 1 in some region 
drives an increase of species 2 in the same region, a clear indication of a condensation type 
of instability. When a varies from zero to ±7r/2 the allowed fluctuations will continuously 
vary from the pure demixing to the pure condensation type. 

For a class of approximations (closures) having the partial direct correlation functions 



vanishing beyond a finite range, it was proven in n| that Q(h) is non-singular for any 
h > 0, so we can limit our search for the zeroes of the minimum eigenvalue to the case 
k — 0. Moreover, since lim^oo det[Q(fc)] = 1 and det[Q(fc)] is a continuous function of h, 
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we must also have det[Q(/c = 0)] non-negative, otherwise det[Q(/c)] would vanish for some 
finite k. We can use this last condition to determine which regions of the phase diagram 
are unstable. We cannot in fact gather this information by just looking at the matrix A, 
which is always positive definite when non-singular. In the following, whenever we omit the 
dependence from the wave- vector k, we shall refer to the case k = 0. 



III. THE BINARY STICKY HARD SPHERE FLUID 



We consider the SHS mixture described in the introduction by the following square- well 



interaction potential between a sphere of species i and one of species j 

+oo < r < a 

-In 

r> R 



14] 



(r) 



1 Rij 

127~2j Rij ~ij 



"13 i 



<Jij <r< R^ , 



(45) 



where (3 = 1/ {k B T) {ks being Boltzmann constant and T the temperature), = (ai+aj)/2 
(<Ti being the diameter of a sphere of species i), R^ — aij denotes the well width, and the 
dimensionless parameter 

- = * = *L > o , (46) 

Tij T T* 

measures the strength of surface adhesiveness or "stickiness" between particles of species 
i and j. In f)4fi)l . r is an unspecified increasing function of T, and we introduced the 
dimensionless quantities = e^/eu and r* = r/e n . The next step which defines the 
SHS model consists in taking the "sticky" limit {Rij} — > {o'ij}. Notice that the logarithm 
in the initial square-well potential (|43|) is chosen so to have a simple expression for the 
Boltzmann factor, which reduces to a combination of an Heaviside step function and a Dirac 
delta function in the sticky limit. 

Within a class of mixed closures for which the partial direct correlation functions Cij(r) 
after the sticky limit vanish beyond [generalized PY (GPY) approximation 2l|], the 
model can be analytically solved for the Baxter factor correlation function 



Qij{ r ) 



\di{r - Uij) 2 + (bi + ai<Tij)(r - a*,-) + K iA , Ly = (a* - aj)/2 <r <a 



, 



elsewhere 



a; = -r 



1 , 120 



A A 2 A 



h = I - a,; ) — 



(47) 



(48) 
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^ n = \ P i<T i ' ^ = J P™ a ™K m , A = 1 - £3 , (49) 

i m 

The Baxter factor matrix Q(/c) first introduced in Eq. ()27j) is related to Baxter factor 
correlation function through 

Qij(k) = 6ij - 27r v /p~p~ %(/V) , (50) 

where %j{k) is the one- dimensional Fourier transform of qij(r). It can be expressed in terms 
of spherical Bessel functions of the zeroth and first order and its explicit expression can be 
found in Eq. (27) of Ref. j^, and will not be reproduced here. 
The symmetric matrix is given by 

K H = Vij > ( 51 ) 

llTij 

where y^ = y%j{<jfj) are the contact values of the partial cavity functions. For this kind of 
system a more natural parameter to use in place of the total density p = J2i Pi is the total 
packing fraction 77 = £3. 

In the modified Mean Spherical Approximation [cjj(r) = fair) when r > Cy, where 



fij(r) = exp [—(3<f>ij(r)] — 1 are the Mayer functions] one can show [2l| that 



Dij = 1 for all i and j , (52) 

In the Percus-Yevick approximation [cij(r) = fii(r)vn(r )} o ne can show that the y~ij have 
to satisfy the following set of coupled quadratic equations 13j 

j/y<Ty = diOij + bi + 2n^2 P k VT~ 1ki(L k i) (53) 

k Tkj 

It is worth stressing that the above expressions are valid for both the mMSA and the 
PY closures, provided that the correct values of y^ are inserted into the matrix given 
in Eq. ()51|) |20l 126 1 . All the results gathered so far in this Section are valid for a generic 
p— component SHS mixture. In the rest of the work we will specialize to two-component 
(p = 2) mixtures. For a binary mixture the determinant of Q(0) can be reduced to the 
following simple expression [30JI 

1 + 2?7 77i Aff + V2X22 



det[Q(0)] 



(l-r/) 2 (l-v) 2 

V1V2 \o/\BT 1 \BT n\BT\ \BT \BT , / \BT\2i / KA \ 
(1 _ 77 )3 ^ A U + A 22 ~ 2A 12 ) - A ll A 22 + ( A 12 ) J > ( 54 ) 
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where 





(55) 



(56) 



BT 

ij 



(1 - ry) A 




(57) 



Our task is the determination of the spinodal line and of the nature of the instability. 
These can be expressed respectively by the reduced temperature r* = f T (paf,xi,(,{^ij}) 
and the angle a = f a (paf, xi, (, {vij}), where ( = oijvi- Sometimes it also proves convenient 
to use another set of independent variables, namely rj,xi,(, {^ij}- 

We anticipate that, while f T will in general depend on the particular chosen closure, f a 
need not mirror this feature. An example is the case studied in Section El where two are 
zero and > 0. Then \j = for i ^ % or j ' ^ j and the spinodal equation 



is sufficient for determining the third A, which turns out to be a function \ij(r), xi, C) in- 
dependent from the particular closure within the class we are considering. Since in each 
matrix element of Q the quantities ijij and r»j appear only in the ratios Ay, it follows that 
the angle a [see Eq. (jJTJ)] will also be independent of the particular closure. 

In the case of a general binary mixture (with two or three non- vanishing e^) we ex- 
pect a dependence of the angle from the closure, even if this point would deserve further 
investigation. 



The PY approximation leads, even in the simple binary case, to the solution of two 
coupled quartic equations. We then start with a simpler task, akin to the one already 
discussed by Chen and Forstmann [5J] for a different potential, of finding the spinodal line 
and angle a predicted by the mMSA and PY for the symmetric binary mixture. In this case 
X\ = X2 = 1/2, 0i = ct 2 = a, and en = €22- By symmetry we must have Cn = £22 and from 
Eq. (|25|) it follows that M is diagonal, the cross term M px being identically zero and 



A_(0) = 



or 



det[Q(0)] = 0. 



(58) 



IV. THE SYMMETRIC BINARY MIXTURE 



A_ =mm{M pp ,M xx } . 



(59) 
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Therefore the symmetric mixture can only have either pure condensation (a = ±7r/2) or 
pure demixing (a = 0) instabilities. 

Moreover for the symmetric mixture we have from Eqs. (|28 )l - (jH()|l 

S PP = xr ' (60) 
lvl pp 

5 » = ' (61) 

S p , = . (62) 

We see then that on a pure condensation instability S pp (0) — > oo or /in(0) + ^12(0) — > 00, 
whereas on a pure demixing instability S xx (0) — > 00 or /iu(0) — ^12(0) — > 00, and each type 
of instability shows a distinct form of long-range behavior in the correlation functions. 



A. Symmetric mixture in the mMSA 

Let us first consider the symmetric mixture within the mMSA. The spinodal line will be 
of pure condensation when M pp = 0, that is 

r' = T ; = (1 + ^l-±, (63) 

whose maximum in the [t*,t]) plane occurs at 77 = <q™ MSA = (\/3 — l)/2 = 0.3660... 
(independently of ^12)- On the other hand the spinodal will be a line of pure demixing when 
M xx = which has as solution 

t* = t: = (1-v 12 )± V . (64) 

Note that the allowed packing fractions are the ones smaller than the close packed packing 
fraction r) = nV2/6 = 0.7404 . . .. 

For the determinant of Q(0) we find from Eq. ()54jl 

det[Q(0)] = ^-^<r Tx) ' + 2 l , (65) 



T 



so that the system is unstable when r lies between the two roots t p and t x , at a given packing 
fraction. 

While the condensation line is always present, the existence of a demixing line depends 
upon the value of z/ 12 , as expected. When V\ 2 > 1 the demixing line r* = r* lies below the 
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?7-axis, and hence the spinodal in the phase diagram (t*,tj) is the curve r* = r* (see Fig. 
12} , with the instability being of pure condensation at all densities. Notice that this would 
be the case for Lorentz-Berthelot mixtures for which we have eu = \f^\\^ri = en, which 
corresponds to u\2 = 1, that is the one-component case. 

When U12 < 1 the two roots t* = r* and r* = r* intercept at a point 31] (see Fig. E} 
having packing fraction 

V = V P x= Q 2UU < 1 , (66) 

3-1/12 

so the instability is of pure condensation for rj < r\ px and of pure demixing for rj > r] px . 



B. Symmetric mixture in the PY 



In the PY approximation we first need to determine the cavity functions at contact. Eqs. 
(J53J) for the binary symmetric mixture can be recast into the following form 



A12T12 



1 + — 



77 1 
r 12 V 2A - -2 VXu 



2A An 



where 



2 + 77 



(67) 
(68) 

(69) 



2{l-r 1 ) 2 ' 

is the HS expression for the cavity functions at contact. Substitution of Eq. (|55jl into Eq. 
(|57} leads to a quartic equation for An. The solution for the cavity functions at contact can 
then be written as 



Til 

V12 

712 



R 



712 



1 + 4(4-^) 



(70) 
(71) 



where R is a solution of the quartic equation. 

In order to find the physically meaningful zeroes of M pp and M xx we proceed as follows. 
First we compute all the four roots Ri, i = a,b,c,d of the quartic equation and hence 
(27ii)i = (y 11) i(r*,r], 2/12), and (3/12)3 = (j7i 2 )i(r*, 77, 1/12) are the cavity functions at contact 
obtained using the root Ri, while (M pp )i and (M xx )i are the diagonal elements of M obtained 
using for the cavity functions at contact (j7n)i and (3/12)1- As it turns out, only two roots 
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Ri will give physically admissible cavity functions at contact. Then we compute the zeroes 
of (M P p)i, denoted as r* = (r*)j(?y, z/ 12 ), and of (M xx )i, denoted as r* = (r*)i(r), Vvi). Then 
physical zeroes are then selected by the requirement that 

\im(y 11 ) i ((T*) i ,r],u l2 )= lim (y 12 )i((T*)i, rj, ^12) = V+ («)<, 17) a = p,x, (72) 

where |/° c is the physical cavity function at contact for the one-component system 

y±(r, 77) = — V : =^ . (73) 



l + i7±V( 1 + l7) 2 -f^ 



Using constraint (J72J) we find that the correct physical solution is R = R b at high 77, the 
only one such that 



lim (y n ) b = lim (y 12 ) 6 = y HS , (74) 

T — >00 T — >00 



\xm(y u ) b = \xm(y 12 ) b = 1 , (75) 

77 — >0 -q— >0 

while at small 77 the solution to use is R = R a such that condition (|72j) is satisfied. As for 
the one-component system there is an interval [0, f] e ] where there are no physical zeroes. For 
the one-component case the spinodal 

„ » I + 477 - 14?7 2 

T T ° C 12(1 - 77) (1 + 277) ' 

exists only if 77 > r] e where r] e = r)f Y = (3v^2 — 4)/2 = 0.1213 . . . and rj^ Y is the PY critical 
packing fraction. For the binary symmetric mixture, numerical results strongly suggests 
the coincidence of rj e with the critical packing fraction (see Figs. EJ and EJ) but we 
have not succeeded in proving it (nor in determining an expression for it). The unphysical 
continuation of the pure condensation spinodal in the range [0, rj e ] is given by the root 
R = R c such that 

lim (yn) c ((r*) c ,T7,i/i 2 ) = lim (£i 2 )c((t*) c , rf, u i2 ) = y -{{T*)c,v) a = p,x . (77) 

viz— ^12—*! 

Notice that this solution would also give, in the same range of 77, an unphysical spinodal of 
pure demixing whenever U12 < 1. 

The zeroes r* = (r*)j and r* = (r*)j are shown in Fig. E]for u± 2 = 2, and for z/12 = 2/3 
(the same conditions as in Figs. El and respectively). As it happened in the mMSA, for 
Ui 2 > 1 there is only a spinodal of pure condensation, while for z/ 12 < 1 a spinodal of pure 
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demixing appears at high 7], as expected on physical grounds. However, unlike the mMSA 
case, the pure demixing and the pure condensation lines do not merge. Also the shapes and 
numerical values of the PY spinodals significantly differ from the mMSA ones. In Fig. El 
we select u%2 slightly above 1 and slightly below 1 in order to check the correct convergence 
towards the one-component case. At v\2 = 1/1.1 the line of pure demixing appears in the 
physically non-accessible region 77 > r] . 



V. TWO PARADIGMATIC SYSTEMS 



The next two mixtures can be regarded as paradigmatic examples of a system where one 
expects to have a predominant condensation or predominant demixing type of thermody- 
namic instability. The thermodynamics of these mixtures have been previously investigated 
by Barboy and Tenne within the PY approximation. In the following we shall extend 
this analysis of the instability type both within mMSA and PY closures. The two systems 
are defined as follows: (System A) ei 2 > 0, en = e 22 = 0; (System B) en > 0, e i2 = e 22 = 0. 
System A corresponds to a fluid where the HS potential acts between like particles and the 
SHS potential between unlike particles, while System B corresponds to HS (species 2) in 
a SHS fluid (species 1). Alternatively, on regarding the large spheres as the "solute" and 
the small spheres as the "solvent" , Systems A and B can be reckoned as a schematic model 
mimicking a "good" and a "poor" solvent, respectively j^ . 

For System A we have A n T = Af 2 T = 0, so Eq. reduces to a a linear equation for y 12 
with the following solution (which corrects Eq. (64) of Barboy and Tenne) 

= irfc ■ (78) 

L ^ 2A ri2 



where 



+ ( 7 q) 

is the HS expression for the contact cavity function. 

For System B we have Af 2 T = \%f = 0, so Eq. (}53|) reduces to a quadratic equation for 
yi\. The only solution which reduces to the HS expression for m — > 00, is (identical to Eq. 
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^ , (80) 



1 + Am + V [ L + A n J 3 2/11 7* 



(57) of Barboy and Tenne) 
Vu = - 

where 

*5' = s + §^- (81) 

is the HS expression for the contact cavity function. The instability lines are again given by 
Eq. (JSED. 

Let Ajj(?7, xi, C) be the solution of the spinodal equation (|58|l for the only non- vanishing 
Ay. As the cavity functions must be positive, the spinodal exists only for those values of 
77, x, ( for which A^j > 0. It may also happen (and it does in the PY case) that the spinodal 
equation 

—(closure)/ j-\ 

Ag we) = ys. — wzi^iij = x _ { ^ Xit j (82) 



upon choosing the correct physical solution for y^f° sure \ does not have any real positive 
solutions for r , at certain values of r], x, (. For these values the spinodal predicted by the 
particular closure has loss of solution and the predicted value for the angle a has clearly no 
physical meaning. 



A. Instabilities for System A 

On setting 




A — 3+ Vl 3+ ^A 3+ ^' (83) 

for System A the solution of Eq. (J58|) within the mMSA approximation is 



2 



r i2 - \bt n n ■ 



while in the PY is 

-ft = v^ MSA - ^ , (85) 

and, in the limit of high dilution while keeping r^ MSA constant, one finds r^ y — > t^ msa , 
as expected in view of the fact that the PY contact cavity functions converge towards the 
mMSA contact cavity functions. 
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In order to exist, the instability line must clearly lie on the r 12 > side of the (r 12 , rj) 
plane. It is easy to see that, while 



j^mMSA 
ar 12 



dr) 

we have that (<7i < a 2 ) 



> for any choice of X\ and ( , (86) 



dr PY 
u/ 12 



> only when — < X\ < „ 2 — ^ . (87) 



dr] 

So in the PY approximation the thermodynamic instability disappears as x\ falls outside 
the range indicated in Eq. (|8*7|) . 

In Figs. El and [7| we depict the mMSA and PY spinodals, respectively, at a given value 
of ( and three different values of X\ for which the PY spinodal does exist. One clearly sees 
that conditions (J86|) and (J87)) result in a large scale difference between the two plots. 

As regards the angle a, we know, from the discussion at the end of Section ITTT1 that the 
angles predicted by the two approximations are the same, and from Fig. |H]we see that the 
kind of instability is mainly of condensation type in accord with what we expected from 
the outset. Two exact limits are worth mentioning. First, the infinite dilute limit 

lim a = arctan ( V^+^M ) , (88) 
n-o \ ^a 2 x 2 /(a 1 x 1 ) - 1 J 

provides an analytical check of the numerical results reported in Fig. |H| Second, when 
<J 2 ^> a i one obtains 



lim a = arctan yx~Jx~2 . (89) 

This result bears an interesting physical interpretation. As the fraction x 2 of large particles 
decreases, the angle a tends to 7r/2, that is to a condensation instability. This is in striking 
contrast with what one would expect for HS on the basis of an entropic depletion mechanism 
(32! which would tend to favor demixing in a system with a small number of large 
spheres. The reason for this can be traced back to the fact that in System A unlike particles 
have attractive interactions, thus preventing smaller particles to slip out from the interstitial 
region between two larger spheres. This interpretation also holds true if one regards system 
A as a good solvent. 
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B. Instabilities for System B 

Denoting 



= (l + 2,)(l-„) 
?7i(l - rj) + 3r?i?72 



for System B the solution of Eq. (|5H|) within the mMSA is 



A 

while in the PY approximation is 



rn MSA = w > ( 91 ) 



.PY _ -HS^mMSA , J7_l ?7l 

jit 
r ll 

for 



r ll — yil r ll T ^2 T mMSA ^ ' ^ JZ / 



^ > 1 (f - l) • (») 

In view of the above constraint, there is an interval rj G [0, rj e ] where no physical spinodal 
exists. We stress that only for the one-component SHS limit (x 2 = 0) one finds that i] e = t] c , 
with rj c being the critical packing fraction, whereas in the more general case, studied here, 
this occurrence is no longer true, as shown in Fig. E3 Once again rfi, as given in Eq. 
()92|) . reduces to t^ msa , in the limit of high dilution, with r^ MSA kept constant. However, 
unlike t™ msa , which is always a concave function of r/ for any choice of x\ and £, r^ y , it 
may display a van der Waals loop (see Fig. fTU|) as a function of r\. The shape of the spinodal 
is strongly dependent on the content of the HS component in the mixture. When x\ < X\ 
[x,\ ~ 0.8681 . . . when ( = 1) the spinodal is a monotonously increasing function of rj, while 
: or x\ > Xi a loop appears. This point has already been emphasized by Barboy and Tenne 

As previously remarked, even in this case both mMSA and PY results for a coincide in 
the respective range of existence. In Fig. ^2 we see that the instability for System B tends 
to pure demixing for ( = 1 and large rj. As ( is increased, one finds the same limit (jH^j) as 
for System A. Once again the osmotic depletion mechanism fails because of the presence of 
stickiness this time among the small particles. As a further support to this interpretation, 
one also finds in the opposite limit 

lim a = arctan I I . (94) 
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In this case, when x\ = r] the instability of the system is of a pure demixing type, so the 
solvent (particles of species 2) is a "poor" one. This is because the smaller particles (species 
2) interact as HS both with larger spheres (species 1) and with each other. Hence, not only 
the depletion mechanism is not opposed in the present case, but, quite on the contrary, is 
favored by the attraction occurring between two big spheres (see Fig. EJ). This results into 
the possibility for the existence of a demixing instability even if the HS binary mixture, 
within the closures considered here, does not have any instability [see Eq. ([98p]. One can 
also show that 

lim a mMSA = arctan ^x 1 /x 2 , (95) 
17 — >o 

where a mMSA is the angle predicted by the mMSA, whose spinodal does not have loss of 
solution at small 77, or, upon using T]i,r]2,y as independent variables, \im m ^ a mMSA = 0. 

Before closing this section, a word of caution should be given on the aforementioned 
interpretations. In order to have a clear and quantitative understanding of the depletion 
mechanism discussed in this section (for both System A and B), the depletion potential, 
that is the effective potential among the large spheres mediated by the presence of the small 
ones, should be computed. Hence, the aforementioned scenarios should only be considered 
as a plausible possibility rather than a definite statement. 



VI. FIVE BINARY MIXTURES TREATED WITH MMSA 



As a final point it is instructive to consider a more general example. To this aim, it 
proves convenient to relate the adhesion strengths eij to the particle sizes Our past 

experience J^J suggests to consider five different cases, obtained setting e^/eo = F^i^i,^) 
for fj, = 1,2, 3, 4, and 5. The functions are selected as follows 



eo 



(<t> 2 /4 Case I , 

ai<jj/crfj Case II , 

(cx^/afj Case III , 

1 Case IV , 



Case V 



(96) 



where (F) = J^Xj-Fj. A critical justification leading to the above choice can be found in 
Ref. [26]. Note that since for all five Cases the are homogeneous functions of order zero in 
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the diameters {cXj}, the corresponding mixtures are invariant under a transformation where 
V — > XV and all <jj — > Act,, with A a scale factor 3^] . 

We have calculated the angle a defined in Eq. fl41j) on the spinodal [Eq. EH] for all the 
Cases listed in ()9(i)l within the mMSA closure. The angle a turns out to be the same for 
Case I and III. The results are shown in Figs. IT2l and ITH1 for x\ = 1/2 and two different 
values of (. We have only considered packing fractions i] < i] m = iry/2/6, where i] m is the 
maximum packing fraction for a "completely demixed" HS mixture (i.e. the packing fraction 
of a mixture where the spheres of species 1 are in a close packed configuration occupying 
a volume V\ and the spheres of species 2 are in a closed packed configuration occupying a 
volume V 2 such that V 2 H V\ = 0). It gives a lower bound to the true maximum packing 
fraction. 

In Cases I and III we have pure condensation as rj — > 0. Case V display a pure condensa- 
tion point at small but non-zero values of 77. In Case II we find a pure demixing point at high 
r), for sufficiently large ( in the same region where in Case IV we have a pure condensation 
point. The packing fraction of pure demixing for Case II can be easily calculated to be 

" = W • (97) 

which turns out to be very close, albeit in general not coincident, with the packing fraction 
at which we find pure condensation in Case IV. 

We remark that (both for mMSA and PY) the presence of an instability curve for the 
SHS model is entirely due to the stickiness, since in the HS limit (r — > 00) we have 

i de « = (i4' (98) 

which is always a positive quantity. Eq. (|98|) can be derived from Eq. ()54|) by noticing 
that the contact values of the partial cavity functions yij must remain finite as r — >• 00. So 
the above statement is actually valid for any closure in which the partial direct correlation 
functions vanish beyond cr^. In particular it is valid for the mMSA and the PY [2^ approxi- 
mations. For other, thermodynamically more consistent closures, the statement is no longer 
rue since phase separation has been observed for highly asymmetric HS binary mixtures 

3- 
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VII. CONCLUSIONS 

In this work we have applied the method devised by Chen and Forstmann [^J to charac- 
terize the kind of thermodynamic instability to a number of carefully selected SHS binary 
systems. The crucial quantity turns out to be the Chen and Forstmann angle a, see Eq. 
(|41|) . on the spinodal: when a is close to the instability is of the pure demixing type, 
whereas a value close to ±7r/2 indicates a pure condensation instability. 

The presence of adhesion between the spheres results in the existence of thermodynamic 
instabilities for the SHS model when treated within closures having the direct correlation 
functions vanishing beyond the hard core ranges, whereas it is known that the HS mixture 
within the same approximations do not show any instability [see Eq. (198)1 ] . 

We have first considered the symmetric binary mixture in the mMSA (see Section II V A)) 
and in the PY approximation (see Section IIVB|) . This latter case was already considered 
by Chen and Forstmann for a different potential. We have found that when en < e 12 the 
instability is of pure condensation along the whole spinodal [see Fig. [2] and Eq. (j63| for 
the mMSA, and Figs. 0] and El for the PY], while when en > ei 2 a pure demixing spinodal 
appears at large packing fractions [see Fig. Eland Eqs. (JH3j) and (JH4*)) for the mMSA, and 
Figs. HJandEJfor the PY], all within their respective limits of validity. This general behavior 
appears to be characteristic of symmetric binary mixtures, in the sense that it is observed 
in systems with pair potentials more "complex" than the SHS potential (hard spheres with 
Yukawa tails [33], square well j^, Lennard- Jones etc.) which do not admit analytic 
solutions. The condensation and demixing lines are found to meet at a point in the mMSA, 
whereas they do not merge within the PY approximation. 

Other two interesting examples can be treated in detail from an analytical point of view 
as discussed in Section We compared the spinodals and the angles a predicted by mMSA 
with those predicted by PY for a binary mixture with e 12 > and en = e 22 = (System A) 
and one with en > and ei 2 = e 22 = (System B). Being the SHS interaction attractive, 
one should expect System A to present mainly condensation instabilities and System B 
mainly demixing instabilities. These choices for the e^ reduce the equations (jHSj) for the 
contact values of the cavity functions in the PY approximation at most to a quadratic 
one, simplifying calculations considerably. We find that the spinodals predicted by the two 
approximations are very different both quantitatively and qualitatively [see Figs. El and 
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and Eqs. flE3) and (jHSJ) for System A, and Fig. [[Hand Eqs. (JHU) and flH2J) for System B]. 
Nonetheless the corresponding angles a do not depend on the closure, when this is chosen 
within the GPY large class containing mMSA and PY as particular cases. In agreement 
with our expectations, we find that the instabilities of System A are predominantly of the 
condensation type (see Fig. |HJ), while the ones of System B of the demixing type when ( ~ 1 
(see Fig. ITT]) . For System B when we have a small number of large spheres of species 1, the 
demixing instability may be favored by both the osmotic depletion mechanism |32 ; ] and the 
stickiness between the large spheres (see Fig. EJ). 

In the more general case, the pair potential depends in general on 3 parameters: the ratio 
of the sphere diameters of the two species, ( = o~2/ 'oi, and two dimensionless parameters 
which measure the relative strength of surface adhesiveness, u 2 2 = 622/611 and = £12/ e n- 
A reduction occurs when the latters are connected to the former through plausible relation- 
ships eij = e§Tij{(T\, a 2 ). Following our previous work [26(, we have considered five possible 
cases [see Section EU and Eq. (l96|)]. We find that four of the five cases exhibit very distinct 
types of instabilities (see Figs. IT2land lT!?|) : Cases I and III have the same angle a, with pure 
condensation at rj — > and predominant demixing for rj > 0; Case V has a pure condensation 
instability point at low packing fractions; Case IV has a pure condensation instability point 
at high packing fractions provided that ( is sufficiently large, whereas Case II has a pure 
demixing instability point under the same conditions. 

It would be desirable to extend the present study in two respects. First it would be 
interesting to consider different, more sophisticated, closures, in view of our results on the 
two examples (denoted as System A and System B) where the angle a is shown to be 
independent of the particular closure within the GPY class, in spite of a large difference 
in the corresponding instability curves. Second, it would be nice to test the analytical 
predictions given in this work against numerical simulations, with a particular attention to 
what concerns the depletion mechanism. We plan to address both issues in a future work. 



APPENDIX A: THERMODYNAMIC RELATIONS FOR THE ELEMENTS OF 

THE M MATRIX 

In this appendix we gather together some well known relationships between thermody- 
namic quantities and the results obtained in the main text. The Ashcroft-Langreth partial 
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structure factors |39| of an homogeneous and isotropic p— component mixture are related to 
the partial total correlation functions as 

Sij(k) = Sij + p y/XiXjhijtk) , (Al) 

where Xi = (iVj) / (N) is the molar fraction of particles of species i and p the total density of 
the mixture. From the normalization condition for the partial pair distribution functions of 
the grand canonical ensemble follows 

The matrix A, defined in Eq. (I13J1 of the text, is related to the structure factors by 

S tJ (k) = [A-'Uk) . (A3) 

We now relate composition fluctuations to thermodynamic quantities. The grand parti- 
tion function is 

oo 

e -W = ^ e ^ P i=i N i^- A (T,v,{m)} ; (A4) 

Ni,...,N p =0 

where A(T, V,{Ni}) is the Helmholtz free energy of a member of the grand canonical en- 
semble with given number of particles of each species, and the chemical potentials {p{\ are 
to be determined from the average number of particles of each species 



(AT.)= N ie ft n+ ZUNm-MT,v,{Ni})] . ( A5 ) 

Ni,...,N p =0 

We immediately find 



) = ~(Ni) , (A6) 



and 



(3 \ dpj J T)V) { K } \9fJ>j / Ty,^} 

= (NiNj) - (NiUNi) = y/xtxjSvmN) , (A7) 

where the index % denotes all species different from i. Since the thermodynamic deriva- 
tives (dNi/dpj)T,v,{fj,j} are the elements of the inverse of the matrix whose elements are 
{dpi/ dNj) T y,{Nj} we can invert the above relation to read 



dN jJT,V,{Nj} ( N )V%i x j 



^4(0) , (A8) 
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where we indicated with S the matrix whose elements are the partial structure factors. 
We now define the partial volumes as 

dV 



(A9) 

T,P,{Ni} 



Since the total volume is an homogeneous function of order one in the extensive variables 
we must have 

v 

Y^N lVl = V, (A10) 



since the Gibbs free energy G = G(T, P, {Aj}) is an homogeneous function of order one in 
the extensive variables we must have 



Y^N tf i t = G, (All) 



i=l 

so in particular the chemical potentials will be homogeneous functions of order zero in the 
variables {Aj}, we can then write = ^(T, P, {{Aj}}) where with the symbol {{A,}} we 
mean that the variables {A;} can appear only as ratios. We also find 

e^A = / dfM_\ + (A12) 

9 n Jt,v, { n j} \ dN J t,p,{{Nj}} V Xt ' 
where xt is the isothermal compressibility 

Notice also that taking the partial derivative of Eq. (jAllJI with respect to Nj at constant 
T, P, and {Aj} we find the following Gibbs-Duhem relation 

tMlfc) =»■ < A14 > 

i=l \ OI "3/ T,P,{{Nj}} 

We want now find thermodynamic relations for the matrix elements M pp , M xx , and M px 
of the binary mixture. We will do the calculation explicitly for M pp and quote the final 
result for the other two elements. So from Eq. (J23j) we find for M pp 

M pp = 2:1(1 - pxicu) +x 2 (l- px 2 c 22 ) - px x x 2 {c X2 + c 2 i) 
= Vp(3 J2 XiXj (^f) 

i,j=l \° ly 3 / T,V,{Nj} 

— / , XiXjViVj 



XT 



(A15) 
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where Xt = P/P * s the isothermal compressibility of the ideal gas, in the second equality 
Eqs. (Unj) and (jA8|) were used, in the third equality we used Eqs. (jA12|) and (jA14|) and in 
the last equality Eq. (jAlOj) . For M px we find 

M px = 6 , (A16) 

XT 

where 

1 fdV\ 

5 = P ( Vl - V2 ) = v ( — ) (A17) 



T,P,N 



and for M, 



XX 



M xx = x 1 x 2 5 2 ^ + x 1 x 2 ^(— Y ) . (A18) 
Xt V \ ox{ 



T,P,N 



2 . (\ r f (d*G 



T,P,N 



, 2 X°t , X°t fd 2 G 

The determinant factorizes 

det(M) = det(A) = [det(Q)] 2 = x x x 2 

Xt v V vXi 
thus yielding Eq. (j2E|) in the main text. 
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LIST OF FIGURES 

Fig. 1 Schematic representation of the two orthonormal vectors z± defined in Eq. (|32|) and 
of the angle a defined in Eq. (14 lj) when a G [0,7r/2]. When a G [— tt/2,0] the angle 
shown in the Fig. corresponds to |a] and 5x to — Sx. 

Fig. 2 Spinodal line (continuous curve) for the symmetric mixture in the mMSA with z/ 12 = 2. 
The kind of instability is of pure condensation along the whole spinodal. 

Fig. 3 Spinodal line (continuous curve) for the symmetric mixture in the mMSA with z/ 12 = 
2/3. In this case the instability is of pure condensation for rj < r] px along r* = r* and 
of pure demixing for 77 > r\ px along r* = r*. 

Fig. 4 Spinodal line for the symmetric mixture in the PY approximation with z/ 12 = 2 in the 
top panel and with vyi = 2/3 in the bottom panel. At z/ 12 = 2 the instability is of 
pure condensation along r* = (r*)j, i = a,b and of pure demixing along r* = (r*)*, 
i = a,b. The zeroes labeled c are unphysical. The gaps between the curves r* = (r*)j 
are numerical artifacts. At z/i 2 = 2/3 there is the appearance of a pure demixing 
spinodal at high 77 which does not cross the pure condensation one. For reference we 
also plot in both panels the spinodal of the one component system r* c [see Eq. (|75j)] 
which is physical only for r\ > r\ c = (3\/2 — 4)/2 = 0.1213 . . .. 

Fig. 5 Same as Fig. 01 with vyi = 1/0.9 in the top panel and z/i 2 = 1/1.1 in the bottom panel. 
In this last case the expected line of pure demixing would start at rj > rjo = 0.7404 . . . 
in the unphysical range of densities. 

Fig. 6 For System A the mMSA spinodal [see Eq. (|84|)] for ( = 2 and three different values 
of X\. 

Fig. 7 For System A the PY spinodal [see Eq. (|85|)] under the same conditions considered in 

Fig. El 

Fig. 8 Behavior of the angle a of Eq. (J41j) predicted by the mMSA and PY for System A when 
x\ = 0.75, C = 2. In this case the PY spinodal has no solutions when rj > 0.03227 . . .. 
In the inset we show the region of rj were the PY spinodal exists. Note that here and 
in the following cos a rather than the angle a itself is depicted for visual convenience. 
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Fig. 9 In System B, when we have a small number of large particles of species 1, the demixing 
instability [see Eq. 1)94)1] should be favored by the osmotic depletion mechanism, since 
the small spheres interact through a HS potential both among themselves and with 
the big spheres. 

Fig. 10 For System B the spinodals predicted by mMSA [thick lines, see Eq. fpTT)) ] and the 
ones predicted by PY [thin lines, see Eq. 1)92)1 ] for ( = 1 at three different values of x\. 
The physically meaningful PY spinodals are those lying above the "existence" lines in 
accord with condition (|93jl . 

Fig. 11 For System B behavior of the angle a of Eq. ()4"T)l predicted by mMSA and PY for 
x\ = 0.91 and £ = 1 in the bottom panel (in this case the PY spinodal has loss of 
solution for rj < rj e ~ 0.1248 . . ..) and £ = 2 in the top panel (in this case the PY 
spinodal has loss of solution for 77 < i] e « 0.1614 . . .). 

Fig. 12 Behavior of the angle a of Eq. ()41)1 for Cases I, II, III, IV, and V when x\ = 1/2 and 
C = 3/2. 

Fig. 13 Behavior of the angle a of Eq. ()41)1 for Cases I, II, III, IV, and V when x\ = 1/2 and 
C = 5. 
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FIG. 1: Schematic representation of the two orthonormal vectors z± defined in Eq. (j.32j) and of 
the angle a defined in Eq. (|41[) when a £ [0, vr/2]. When a £ [—ir/2, 0] the angle shown in the Fig. 
corresponds to \a\ and 5x to — 5x. 
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FIG. 2: Spinodal line (continuous curve) for the symmetric mixture in the mMSA with 1/12 = 2. 
The kind of instability is of pure condensation along the whole spinodal. 




FIG. 3: Spinodal line (continuous curve) for the symmetric mixture in the mMSA with v\2 = 2/3. 
In this case the instability is of pure condensation for rj < r) px along r* = r* and of pure demixing 
for 77 > f] px along r* = r*. 
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FIG. 4: Spinodal line for the symmetric mixture in the PY approximation with 1^12 = 2 in the top 
panel and with U12 = 2/3 in the bottom panel. At v±2 = 2 the instability is of pure condensation 
along r* = (r*)j, i = a,b and of pure demixing along t* = (r*)j, i = a,b. The zeroes labeled c are 
unphysical. The gaps between the curves r* = (r*), are numerical artifacts. At v\2 = 2/3 there is 
the appearance of a pure demixing spinodal at high rj which does not cross the pure condensation 
one. For reference we also plot in both panels the spinodal of the one component system r* c [see 
Eq. COIl] which is physical only for 77 > 7] c = (3^2 - 4)/2 = 0.1213 . . .. 




FIG. 5: Same as Fig. 0Jwith V\% = 1/0.9 in the top panel and v\2 = 1/1.1 in the bottom panel. 
In this last case the expected line of pure demixing would start at rj > r]Q = 0.7404. . . in the 
unphysical range of densities. 
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FIG. 6: For System A the mMSA spinodal [see Eq. (|84j) ] for £ = 2 and three different values of x\. 
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FIG. 7: For System A the PY spinodal [see Eq. (|H51 )] under the same conditions considered in 
Fig. El 
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FIG. 8: Behavior of the angle a of Eq. (|41|) predicted by the mMSA and PY for System A when 
x\ = 0.75, C = 2. In this case the PY spinodal has no solutions when r) > 0.03227 .... In the inset 
we show the region of rj were the PY spinodal exists. Note that here and in the following cos a 
rather than the angle a itself is depicted for visual convenience. 





FIG. 9: In System B, when we have a small number of large particles of species 1, the demixing 
instability [see Eq. (|94[) ] should be favored by the osmotic depletion mechanism, since the small 
spheres interact through a HS potential both among themselves and with the big spheres. 
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FIG. 10: For System B the spinodals predicted by mMSA [thick lines, see Eq. (|91|l] and the ones 
predicted by PY [thin lines, see Eq. I|92jl] for £ = 1 at three different values of x\. The physically 
meaningful PY spinodals are those lying above the "existence" lines in accord with condition (|93|) . 
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FIG. 11: For System B behavior of the angle a of Eq. (|41[) predicted by mMSA and PY for 
x\ = 0.91 and £ = 1 in the bottom panel (in this case the PY spinodal has loss of solution for 
i] < i] e ~ 0.1248 . . ..) and £ = 2 in the top panel (in this case the PY spinodal has loss of solution 
for r] < 7] e f« 0.1614. . .). 
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FIG. 13: Behavior of the angle a of Eq. (|41|) for Cases I, II, III, IV, and V when x\ = 1/2 and 
C = 5. 



